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Abstract. With the invention of new far- infrared (FIR) and radio mm and sub-mm 
instruments (DESIR on SMESE satellite, ESO-ALMA), there is a growing interest in 
observations and analysis of solar flares in this so far unexplored wavelength region. 
Two principal radiation mechanisms play a role: the synchrotron emission due to ac- 
celerated particle beams moving in the magnetic field and the thermal emission due to 
the energy deposit in the lower atmospheric layers. In this contribution we explore the 
time-dependent effects of beams on thermal FIR and radio continua. We show how and 
where these continua are formed in the presence of time dependent beam heating and 
non-thermal excitation/ionisation of the chromospheric hydrogen plasma. 
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Introduction 



New ground-based and space-borne instruments are being designed to ob- 
serve the Sun in so-far unexplored wavelength windows or with unprece- 
dented spatial and/or temporal resolution. Here we are motivated by the 
planned observations of solar flares in far infrared (FIR) region and in sub- 
millimetre/millimetre (SMM) radio wavelengths. A specific FIR window 
where the radiation is not observable from the Earth is between 35 and 
250 /im and this is expected to be explored by the DESIR te lescope on 



board the SMESE satellite (SMall Explorer for Solar Eruptions - I Vial et al 



( 200 7])) . At SMM wav elen gths, observations were already performed by 
Kaufmann et al.l feOQOh anc ULuthi et~ail J20041 ). However, new data are ex- 



pected soon from ALMA (Atacama Large Millimeter Array) radio-interferomete 
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In both FIR and SMM domains a superposition of two radiation components 
is expected: the high-frequency part of the microwave spectrum from ultra- 
relativistic electrons and/or positrons (non-thermal synchrotron emission) 
an d the low-frequen cy part of the thermal flare continuum. As suggested 



m 



Vial et al 



20071 ). these emissions provide, respectively, a unique diag- 
nostic of the non-thermal particles (beams) and of the thermal response 
of the lower flare atmosphere to the energy deposit. In FIR domain, the 
relativ e importan ce of both process e s was esti mated by Qhki and Hudson! 
J197RI ). Recently, lLonkitcheva et all hoo4 ) and iHeinzel and Avrettl (l2009h 
computed the microwave thermal continua formed in dynamic and static 
semi-empirical atmospheres, respectively 

In this study, we concentrate on formation of such thermal emission in 
the 35 /um - 1 cm wavelength range originating in a highly dynamical flare 
atmosphere heated by short-duration electron beam pulses. Section 2 de- 
scribes the mechanisms responsible for the thermal emission, characteristics 
of the emission corresponding to simulations of the beam heated atmosphere 
are presented in Section 3. Section 4 summarises our results. 



2. Mechanisms of the FIR and SMM thermal continuum 

formation 

Depending on the wavelength, various atm ospheric depths can contri bute 
to the emergent FIR and SMM intensity (Irleinzel and Avrettl . l2009l l. In 
the chromosphere, the dominant source of opacity i s the hydrogen free-free 
conti nuum. The absorption coefficient is given as (IRybicki and Lightmanl . 
197$ ) 

Kj/(H) = 3.7 x W 8 T- 1/2 n c n p v- 3 g s (1) 

where n e and n p are the electron and proton densities, respectively, T is the 
temperature, and gs sa 1 is the Gaunt factor. At lower temperatures, around 
the temperature minimum region, H~ free -free opacity plays an important 
role. The H~ absorption is (iKnnicd . ll97oh 



k„(H- 



n e n H 



v 



(1.3727 x 10" i& + (4.3748 x 10~ iU -2.5993 x 10" 7 /T) jv) (2) 



where «h is the neutral hydrogen density. Total absorption corrected for 
stimulated emission is then 



[k„(H) + k„(H-)] (1-e 



-hu/kT\ 



(3) 



2 
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where h and k are the Planck and Boltzmann constants, respectively. The 
source function for these free- free processes is the Planck function B v . How- 
ever, the absorption has to be evaluated using number densities coming 
from non-LTE calculations. Finally, emergent intensity I u can be obtained 

as 

I v = I r] u e~ Tv dz r\ v = k v B v &t p = -K u dz , (4) 



where z is the geometrical height. 

3. Flare models 

We have evaluated the FIR and SMM continuum emission for non-LTE 
radiative hydrodynamic models which describe the atmospheric response 
to the time dependent electron beam heating. In these simulations elec- 
tron beams of a varying beam flux and power-law spectra are injected into 
a loop which correspo nds to an initially hydrostatic VAL C atmosphere 



(jVernazza et all 119811 ). Then, the time evolution of the beam particles and 
plasma properties are followed. We also consider the so-called non-thermal 
collisional rates, C nt , which account for the changes of the atomic level 
population s due to the collisio n s with the beams . For m ore details on the 
models see IVaradv et all fcOQSl ): iKasparova et al.l (j2009l ). 



3.1. Time variation 

Our simulations show that the continuum emission is well correlated with 
the beam flux on the time scale of the beam flux variation, e.g. on a sub- 
second time scale. Longer wavelengths (1 cm) exhibit larger relative inten- 
sity increase than shorter wavelengths (35 //m) when compared to the the 
preflare (VAL C) intensities - see Figure [U Furthermore, emission in the 
A < 0.2 mm range depends also on the beam parameters, i.e. the power- 
law index and the beam flux. Generally, maxima of the emission lag behind 
the beam flux maxima, the time lag lies within the 0.1 - 0.3 s interval for 
present simulations. 

3.2. Influence of non-thermal collisional rates 

Although the non-thermal rates can significantly affect the hydrogen level 



populations and electron densities in some atmospheric layers (jKasparova et al.l . 
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Figure 1: Time evolution of relative continuum intensities scaled to the intensity at t = s, 
i. e. of the VAL C atmosphere, for an electron beam of a power-law index <5 = 3 and a 
maximum beam flux F max = 6 x 10 10 erg cm" 2 s" 1 . The dashed line indicates the beam 
flux time variation. 



20051 ). their influence on FIR and SMM continuum emission is almost neg- 
ligible. C Ilt cause only moderate and temporary increase of the continuum 
intensity shortly after the beam injection into the VAL C atmosphere. Ad- 
ditionally, influence of C nt for A < 0.2 mm depends on the beam parameters 
and it is modulated by the beam flux time variation - see Figure [2l 



3.3. Formation depths 

Atmospheric layers which contribute the most to the outgoing intensity, i.e. 
the formation depths, can be understood in terms of the so-called contribu- 
tion function CF 

CF = r] u e- Tv I u = [ CFdz, (5) 
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Figure 2: Influence of the non-thermal collisional rates C nt . The plot shows the ratio of 
intensities with and without C nt . The dashed line indicates the beam flux time variation. 
The beam parameters are the same as in Fig. [T] Note the reversed wavelength scale with 
regard to Figure [T] 



compare with Eq. (jlj). In the initial atmosphere, emission at wavelengths 
shorter than A ~ 1 mm comes from deep, photospheric layers, whereas the 
continuum at other wavelengths is formed at the upper part of the atmo- 
sphere. However, as the atmosphere is heated by the beams, the emission 
is formed mainly at the chromospheric heights (above z ~ 1000 km). The 
emission which comes from these layers can be influenced also by C nt since 
the beams deposit their energy there. On the other hand, the emission in the 
35 - 200 /im range partially originates from the photosphere - see Figure 03 
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Figure 3: An example of the time evolution of contribution function CF (Eq. [5} for an 
electron beam of 5 = 3 and F max = 4.5 x 10 10 erg cm" 2 s" 1 . Left: CF for the VAL 
C atmosphere, i.e. at t = s. Right: CF at t = 1.5 s corresponding to the model with 
included C nt . Colour scale going from white to black represents the range of log CF values 
(from minimum to maximum). Black solid curve indicates the total energy deposit of the 
beam, dashed curve is for the energy deposit on hydrogen En- C nt are proportional to 
Eh, therefore the dashed line marks the layers which are affected by C nt . Gray lines 
display the levels of optical depth r = 0.1, 1, 10. 



4. Conclusions 

Using the results of non-LTE radiative hydro dynamic simulations we have 
found that the rapidly varying beam flux can manifest itself in the thermal 
FIR and SMM continuum intensities. We show that the continuum inten- 
sities do vary on the beam flux variation time scales. However, the direct 
influence of the beam electrons via the non-thermal collision is moderate 
only. The emission is mainly affected by the temperature variations which 
result from the time dependent heating by the beams. Concerning the for- 
mation of the emission, in the heated atmosphere or due to the non-thermal 
collisions the continua are formed in the upper atmosphere. The intensities 
in the 35 - 200 |im partially originate also from the photosphere. 

In this work we did not consider effects related to the situation when the 
frequency of the emitted electromagnetic waves is close to the local plasma 
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frequency. In the case of the solar atmosphere plasma, these effects must be 
taken into account at the radio frequency range. We will address this issue 
for the FIR and SMM emission in a future study. 
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